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Abstract 

We calculate direct gaugino pair production at hadron colliders at next-to-leading order of per- 
turbative QCD, resumming simultaneously large logarithms in the small transverse-momentum and 
threshold regions to next-to-leading logarithmic accuracy. Numerical predictions are presented for 
transverse momentum and invariant mass spectra as well as for total cross sections and compared 
to results obtained at fixed order and with pure transverse-momentum and threshold resummation. 
We find that our new results are in general in good agreement with the previous ones, but often 
even more precise. 
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I. INTRODUCTION 



Weak-scale supersymmetry (SUSY), and in particular the Minimal Supersymmetric Stan- 
dard Model (MSSM), is a theoretically and phenomenologically well- motivated extension of 
the Standard Model (SM) of particle physics [1]. Consequently, the experimental search for 
the spin partners predicted by the MSSM for each of the SM particles is one of the defining 
tasks at current high-energy colliders such as the p»_Tevatron collider at Fermilab |2| and 
the pp Large Hadron Collider (LHC) at CERN 0, |4|. Of particular interest are the neu- 
tral and charged fermionic partners of the electroweak gauge and Higgs bosons, which mix 
due to their equal quantum numbers into four neutralino (x°, i = 1, ... ,4) and chargino 
i — 1)2) mass eigenstates, since these participate virtually always in SUSY collider 
signatures and may also have important implications for dark matter and cosmology. Their 
decays into leptons and missing transverse energy, carried away by the lightest SUSY particle 
(LSP, often the x?), are easily identifiable at hadron colliders. The lighter gaugino/higgsino 
mass eigenstates are accessible not only at the LHC with center-of-mass energies y/S of 7 to 
14 TeV js, 4]], but also at Run II of the Tevatron (\/S =1.96 TeV), where the production 
of X1X2 pairs decaying into trilepton final states is one of the gold-plated SUSY discovery 
channels 

For an efficient suppression of the SM background from vector-boson and top-quark pro- 
duction and a precise determination of the underlying SUSY-breaking model and masses, an 
accurate theoretical calculation of the signal (and background) cross section is imperative. 
As the LSP escapes undetected, the key distribution for SUSY discovery and measurements 
is the missing transverse-energy ($r) spectrum, which is typically restricted in the experi- 
mental analyses by a cut of 20 GeV at the Tevatron and 30 GeV at the LHC. While the SUSY 
particle pair is produced with zero transverse momentum (p?) in the Born approximation, 
the possible radiation of gluons from the quark- ant iquark initial state or the splitting of glu- 
ons into quark-antiquark pairs at 0(a s ) in the strong coupling constant induces transverse 
momenta extending to quite substantial values and must therefore be taken into account. 
In addition, the perturbative calculation diverges at small pt, indicating the need for a 
resummation of soft-gluon radiation to all orders. Only after a consistent matching of the 
perturbative and resummed calculations an accurate description of the (missing) transverse 
energy spectrum and precise measurements of the SUSY particle masses can be achieved. 
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Furthermore, when the SUSY particle pair with invariant mass M is produced close to the 
production threshold at the partonic center-of-mass energy s, soft gluon emission leads again 
to potentially large (logarithmic) terms, which must be resummed to all orders in order to 
obtain a reliable cross section. The production of SUSY particles at hadron colliders has 
been studied at leading order (LO) of perturbative QCD since the 1980s |5fl. More recently, 
previously neglected electroweak contributions jfj, polarization effects [j], and the violation 



of flavor 



and CP symmetry jsj have been considered at tins order. Next-to-lea.ding order 



(NLO) corrections have been computed within QCD since the late 1990s lOj and recently 



also within the electroweak theory 111. Resummation at the next-to- leading logarithmic 

nn 

(NLL) level has been achieved in the small-p^ region for sleptons and gaug inos jl 21 . |13| and 



14 



in the threshold region for sleptons, gauginos, squarks and gluinos 

Since the dynamical origin of the enhanced contributions is the same both in transverse- 
momentum and threshold resummations, i.e. the soft-gluon emission by the initial state, it 
would be desirable to have a formalism capable to handle at the same time the soft-gluon 
contributions in both the delicate kinematical regions, pt <C M and M 2 ~ s. This joint 



resummation formalism has been developed over the last twelve years 16]. The exponen 



tiation of the singular terms in the Mellin (N) and impact-parameter (b) space has been 
proven, and a consistent method to perform the inverse transforms, avoiding the Landau 



po 



e and the singularities of the parton distribution functions (PDFs), has been introduced 
17J • Applications to prompt-photon, electroweak and Higgs boson, heavy-quark and slepton 
pair production at hadron colliders have exhibited substantial effects of joint resummation 
on the differential and total cross sections 

In this paper, we present the first calculation of joint resummation of soft gluon effects in 
the small transverse momentum and threshold regions for gaugino/higgsino hadroproduction 
at the NLL level, using the formalism described above. As in our previous calculations 
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151 ]. we include not only the QCD, but also the SUSY-QCD virtual loop contributions with 



internal squark mixing in the hard coefficient function of the resummed cross section, which 
therefore reproduces, when expanded and integrated over p?, the correct NLO SUSY-QCD 
cross section in the threshold region. For the Tevatron, we consider not only the production 
of xfx®i but also of X2X2 an d X1X1 pairs. In particular the latter can have significantly 
larger cross sections than trilepton production due to the s-channel exchange of massless 
photons. For the LHC, we concentrate on predictions for its current center-of-mass energy 
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of \/~S =7 TeV, where threshold effects and direct gaugino pair production (as opposed to 
the production from squark and gluino cascade decays) will be more important. However, 
we will also show cross sections for the production of light (xf, x°) an d heavy (xf , X34) 
gaugino combinations at the LHC design energy of y/S =14 TeV for completeness. 

The remainder of this paper is organized as follows: In Sec. [Ill we briefly review the 
joint resummation formalism in Mellin and impact parameter space, giving the explicit 
form of the resummed logarithms at NLL order as well as their matching to the fixed 
order calculation and the prescriptions employed for the inverse integral transforms. Sec. II III 
contains numerical results for various gaugino pair production cross sections at the Tevatron 
and at the LHC and compares the results obtained in joint resummation to those obtained 
with the pure px- and threshold resummation formalisms, respectively. We summarize our 
results in Sec. HV] 



II. JOINT RESUMMATION FORMALISM 

Thanks to the QCD factorization theorem, the double differential cross section for the 
hadronic production of a final state with fixed invariant mass M and transverse momentum 
Pt 

d 2 cr 

M 2 — -^^"(r) = ^2 dx a dx b dz[x a f a/A (x a ,^)][xJ b/B (x b ,fip)][zda ab (z,M 2 ,pl,fi 2 F )} 

T ab 

x 8(t - x a x b z) (1) 

can be obtained by convolving the partonic cross section da ab with the universal densities 
fa,b/A,B of the partons a, b, carrying the momentum fractions x a ,b of the colliding hadrons 
A,B, at the factorization scale y,p. The application of a Mellin transform 



F(N) 



fdyy N ^F{y) (2) 
Jo 



to the quantities F G {(Jab, a ab , f a/A , f b/B ] with y G {r = M 2 /S, z = M 2 /s, x a , x b } 
allows to express the hadronic cross section in moment space as a simple product, 

M2 dM^dp 2 T {N ~ l) = E^^'^^^^^l^M 2 ,^,^). (3) 

Furthermore, the application of a Fourier transform to the partonic cross section a ab allows 
to correctly take into account transverse-momentum conservation, so that in moment (N) 





FIG. 1: Tree-level Feynman diagrams for the production of gaugino pairs, 
and impact parameter (b) space it can be written as 



a ab {N, M 2 ,p 2 T ,^ 2 F ) 



db-J (bp T )a ab (N, M 2 ,b 2 ,n 2 F ). 



Here, Joiy) denotes the O th -order Bessel function and 

oo 

a ab (N,M 2 ,b 2 ,fi 2 F ) = ^a?(rfJa%\N,ltf, 



(4) 



(5) 



n=0 



is usually expanded perturbatively in the strong coupling constant a s (fi 2 ) = a s (/i 2 )/(27r) at 
the renormalization scale /ip. For simplicity, we identify in the following the factorization 
and renormalization scales, i.e. fip = = ^- 

In the Born approximation, the hadropro duct ion of neutralinos xt an d charginos xt is 
induced by quarks q and antiquarks q' in the initial (anti-)protons and is mediated by s- 
channel electroweak gauge-boson and t- and w-channel squark exchanges (see Fig. [p. Its 
partonic cross section a^, can be expressed in terms of the gaugino and squark masses 
m-o,± and m^, the masses of the electroweak gauge bosons, and generalized charges {jj. At 



0(a s ), virtual loop and real parton emission corrections must be taken into account 10]. 
The latter induce not only a deviation of the partonic center-of-mass energy s from the 
squared invariant mass M 2 of the gaugino pair, but also non-zero transverse momenta px, 
that extend typically to values of the order of the gaugino mass. 

Close to the partonic production threshold, where z = M 2 /s — > 1 or iV — > oo, the 
convergence of the perturbative expansion is spoiled due to soft gluon radiation, which 
induces large logarithms 

'ln m (l-z) 



with m < 2n - 1 and N = Ne lE 



n i„m+l 



->■ a™ In 



N + . . . 



(6) 



151 ] . Similarly, in the small-p^ (or large-6) region, where 



the bulk of the events is produced, the convergence of the perturbative expansion is again 
spoiled by soft gluon radiation, which induces large logarithms 

1 , „ M 2 



^ln r 
Pt 



-5- ) -> a^\n m+l b z + ... 
Pt 



(7) 



with m < 2n — 1 and b = bMe 7E 12 13 1. The crucial observation, first made by Li 161 and 

nn 

then further developped by Laenen, Sterman, and Vogelsang [17|, |18j is that the common 
kinematic origin of these divergences allows for a joint resummation of the large logarithms 
in the partonic cross section by choosing a function 

N 



(8) 



1 + T]b/N 

which interpolates between N in the threshold region, N 3> 6, and b in the small-pr region, 
b ^> N. Its exact form is constrained by the requirement that the leading and next-to- 
leading logarithms in b and N are correctly reproduced in the limits b — > oo and N — > oo, 
respectively. The choice of Eq. (JBJ) with r\ > (we use 77 = 1) avoids the introduction of 
sizeable subleading terms into perturbative expansions of the resummed cross section at a 
given order in a s , which are not present in fixed-order calculations 18]. Up to NLL order, 
the resummed cross section can then be written in the form 

a ab (N,M 2 ,b 2 ,^) = 4 1 i(^^/x a ,M 2 )^S(^^/x 2 ,A* a )C'^(iV ) o.(^/x 2 )) 

a',a",V,V 

x Cy lv (N,a s (M 2 /x 2 ))H aliyi (M\^)eMGa^(M\N;b\^)}, (9) 

which is very similar to the one of the py-resummed cross section. The operators E a i a ^b allow 
to evolve the PDFs f a ,b/A,B{N, n 2 ) from the scale \j? to the scale M 2 / y 2 . They satisfy, like the 
PDFs themselves, the Altarelli-Parisi equations 19J] and can be written to one-loop order, 
and in the singlet /non-singlet basis in closed exponential form 20j. When the finite 
part of the renormalized virtual one- loop contribution Ao, which in our calculation contains 
the full SUSY-QCD corrections with internal squark mixing [la ], is entirely absorbed into 
the hard function 



H ab (M 2 ,fi 2 ) = af b \M 2 ,f, 2 )[l + a s Ao}+0(a 2 s ), 



(10) 



the coefficients C ab become process-independent and can be written up to one-loop order as 

^(o) 

'ah 



C:°J = 5 ab and 



C$(N) = C a -S ab -P^(N), 



7l~ 



(11) 



where C ab = X^^Lo a s ^1* j ^ ne QCD color factors are C q = Cp = 4/3 and C g = Ca = 3, 
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respectively, and 

p(iWjV) = ~° F pW> e (N) = ~ 2Tr (12) 

11 1 > N(N + 1) ' is v > (JV + l)(JV + 2)' 1 ' 

P^(N) = j^L , p(J), W = (13) 

represent the C?(a s , e) terms in the Altarelli-Parisi splitting functions in Mellin space with 
Tr = 1/2. The important feature of Eq. (Q is that the hard function H a i,(M 2 , /z 2 ) , contrary 
to the perturbative cross section a a b(N, M 2 ,b 2 , fi 2 ) in Eq. (jSJ), no longer contains large log- 
arithms in either N or b, since these have all been resummed in the coefficient functions C, 
the evolution operators E, and the Sudakov exponent G. The latter can be expanded as 

G a b(M 2 , N,b 2 , fj 2 ) = Lg^(\)+g%(\MN,M 2 /^ 2 )+ ••• (14) 

where L = ln% = A/(as/3 ) encodes the large logarithms in both N and 6, /3 = — 
2NfTfj/3 is the one-loop coefficient of the QCD beta-function, and Nf is the number of 
active quark flavors. The first term in this expansion 

Lg^W = ^r(4? } + [2A + ln(l - 2A)] (15) 
collects the leading logarithmic (LL) contributions, while the second term 
2ft, S <?(A,lniV,MV) = (AW + 4») [2A 1 ~ ff^ S + ln(l - 2A)] In 

+ (4" + 4")|[(2A + Ml - 2A)) i-^fti°^ + i lD 2 (1 _ 2A)] 

- (4 2 '+4 2 ')^[2A 1 "^° A ln ^ + ln(l-2A)] 

+ (-2 7 <" - 2 7 « + C») ln(l - 2A) (16) 

with 0i = {11C\ — 5CaNj — 3CpNf)/6 collects the NLL contributions. These terms re- 
produce those appearing in the transverse-momentum resummation formalism [k^, when 



b > N, x -> b and (1 - 2a s /3 In N 
threshold resummation formalism 



JJA — 2 A) — > 1/(1 — 2A), as well as those appearing in the 



15|, when N > b, x -> N and (l-2a s /3 ln iV)/(l-2A) -> 



1, respectively. For transverse-momentum resummation, the logarithm (In b) had to be mod- 
ified (to In a/1 + b 2 ) in order to suppress unphysical resummation contributions at large px 
and small b. This is not necessary within the joint resummation formalism, as the small-6 
singularity is regularized through the function x(iV,6). The coefficients 



A® = 2C a , = 2C a 



and D§ = (17) 
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are well-known from both threshold and transverse-momentum resummation 2l(, while the 
anomalous dimensions 

3C F 



2 ' 7f? 

of the quark and gluon fields have been introduced in order to remove the corresponding 
NLL terms from the one- loop approximation of the diagonal evolution operators E^J . 

While the large logarithms must clearly be resummed close to the production threshold, 
when z — > 1 and N — > oo, and/or at small values of p F — > 0, when b — > oo, they account only 
partially for the full perturbative cross section away from these regions. In order to obtain 
a valid cross section at all values of z and pt, the fixed-order (f.o.) and the resummed (res.) 
calculations must be combined consistently by subtracting from their sum the perturbatively 
expanded (exp.) resummed component, 



(i) 



0o 



(18) 



^(res-) , _(f-0.) _ ^-( ex P-) 
°ab °» 



— a 



ab 



(19) 



The latter is easily obtained by expanding Eq. (jlj) to the desired accuracy. At 0(a s ), one 
finds 



2 ^- ln ^T ) E W ) (M 2 ,^ 2 )^ 1) (iV) + ^)(iV)^(M 2 ,^)] 



+ a s J2 [H^{M\v?)Ci\N) + C£){N)H%\M\n 

c 

- a s H^{M\^)[j\A^ + A$>) - 2J( 7 « + 7 ; i} 
where the O(a s ,e ) terms of the Altarelli-Parisi splitting functions are 

N 





= C F 




= C F 




= T R 


P£\N) 


= 2C, 



1 X ^ 1 



2 N(N + 1) 

2 + N + N 2 
N(N 2 - 1) 
2 + iV + N 2 
_N(N + l)(N + 2) 

1 



1 



_N(N-1) (N + l)(N + 2) f^k 



N ^ 



+ /3o 



(20) 

(21) 

(22) 
(23) 

(24) 



and the full dependence on the transverse momentum p F is embodied in the Bessel integral 

b 

J = I db-J (bp T )\n X (N,b). (25) 



o 
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The form of this integral has been chosen to match the one encountered in pT-resummation. 
Unfortunately, it then induces 7r 2 -terms which differ slightly from those encountered in 
threshold resummation jl8 |. While these terms are formally beyond NLL order, one should 
nevertheless expect slightly better numerical agreement with the p^-resummed calculation 
than with the threshold-resummed calculation. 

After the resummation of the partonic cross section has been performed in N- and b- 
space, we have to multiply the resummed cross section and its perturbative expansion with 
the moments of the PDFs f a /A{N, /i 2 ) and transform the hadronic cross section obtained in 
this way back to the physical z- and p^-spaces. The moments of the PDFs are obtained 
through a numerical fit to the publicly available PDF parameterizations in x-space. For 
the inverse integral transforms, special attention has to be paid to the singularities in the 
resummed exponents, i.e. when A = 1/2 in Eqs. f|T5|) and f fT6|) . They are related to the 
presence of the Landau pole in the perturbative running of a s (/z 2 ), and prescriptions for both 
the Mellin and Fourier inverse transforms are needed. For the Fourier inverse transform of 



Eq. (p[J), we follow Ref. 17j] and deform the integration contour of the b- integral in the 



complex plane by defining two integration branches 

b = (cos ip ± i sin (p)t, £e[0,oo[, (26) 

where (p has to be chosen in the range ]0, 7r/2[. The Bessel function Jo(y) in Eq. fll]) is then 
replaced by the sum of the two auxiliary functions 



h(y,v) = ~ I d6e~^\ (27) 

^ J-iVTT 

-IVTT 



h 2 (y,v) = ~ I dOe-'y^, (28) 

which are finite for any value of y. Their sum is independent of v and is always equal 
to Jo(y)- Since the two functions distinguish positive and negative phases in the complex 
6-plane, they can be associated with only one of the two branches. For the inverse Mellin 
transform 

F(y) = [ ^y- N F(N) (29) 

Jc N 2m 

we choose an integration contour Cn according to the principal value procedure proposed in 



Ref. [22] and the minimal prescription proposed in Ref. 23[ and define again two branches 

C N : N = C + ye ±i *, ye[0,oo[. (30) 



The parameter C must be chosen in such a way that the poles in the Mellin moments 
of the parton densities, which are related to the small- x (Regge) singularity f a /A{ x i lA) (X 
x a (l — x)P with a < 0, lie to the left and the Landau pole to the right of the integration 
contour, respectively. While formally the angle <ft can be chosen in the range [7r/2,7r[, it is 
advantageous to take <fi > it/2 to improve the convergence of the inverse Mellin transform. 



III. NUMERICAL RESULTS 



We now turn to our numerical analysis of joint resummation effects on the production 
of various gaugino pairs at the Tevatron pp-collider (\/S = 1.96 TeV) and the LHC pp- 
collider (VS = 7 - 14 TeV). For the masses and widths of the electroweak gauge bosons, 



we use the values of m z = 91.1876 GeV and m w = 80.403 GeV [24|. The CKM-matrix 
is assumed to be diagonal, and the top quark mass is taken to be 173.1 GeV 25j]. The 
strong coupling constant is evaluated in the one-loop and two-loop approximation for LO 
and NLO/NLL+NLO results, respectively, with a value of A^T 5 corresponding to the em- 

1V1 o 

ployed LO (CTEQ6L1) and NLO (CTEQ6.6M) parton densities 26j. For the resummed and 
expanded contributions, the latter have been transformed numerically to Mellin iV-space. 
When we present spectra in the invariant mass M of the gaugino pair, we identify the un- 
physical scales fip — — A* with M, whereas for transverse momentum distributions and 
total cross sections we identify them with the average mass of the two produced gauginos. 
The remaining theoretical uncertainty is estimated by varying the common scale [i about 
these central values by a factor of two up and down. The running electroweak couplings 
as well as the physical masses of the SUSY particles and their mixing angles are computed 



with the computer program SPheno 2.2.3 



271 ]. which includes a consistent calculation of 



the Higgs boson masses and all one-loop and the dominant two-loop radiative corrections in 
the renormalization group equations linking the restricted set of SUSY-breaking parameters 
at the gauge coupling unification scale to the complete set of observable SUSY masses and 
mixing angles at the electroweak scale. 

For the Tevatron, we choose the low-mass point LMO (SU4) with universal fermion mass 
mi/2 — 160 GeV, scalar mass m = 200 GeV, trilinear coupling A = —400 GeV, bilinear 
Tiggs mass parameter fi > 0, and ratio of Higgs vacuum expectation values tan/3 = 10 
ft nas been been defined by the CMS (ATLAS) collaboration with the objective of 
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high cross sections and thus early discovery at the LHC, as the resulting gaugino, squark 
and slepton masses m^o = m-± = 113 GeV, m^o = 61 GeV, m„ ~ 420 GeV, and mr o± 220 

^\,2 ^ 1 ^1 

GeV lie just beyond the current Tevatron limits. In this scenario, the lightest chargino and 
second-lightest neutralino decay with 35% and 15% probability through virtual sleptons to 
the LSP and one and two charged leptons, respectively 2o| . 

For the LHC, we choose the widely used minimal supergravity (mSUGRA) point SPSla' 



30[ as the benchmark for our numerical studies. This point has the same intermediate 



value of tan/3 = 10 and /x > (favored by the rare decay b — > 57 and the measured 
anomalous magnetic moment of the muon), a still relatively small gaugino mass parameter 
of mi/2 — 250 GeV, and a slightly lower scalar mass parameter m = 70 GeV and trilinear 



coupling A = —300 GeV than the original point SPSla [3l| in order to render it compatible 
with low-energy precision data, high-energy mass bounds, and the observed cold dark matter 
relic density. It is also similar to the post-WMAP point B' (m = 60 GeV and A = 0) 
32]. which has been adopted by the CMS collaboration as their low-mass point LM1 jj]. 
In the SPSla' scenario, the \i is the LSP with a mass of 98 GeV, the gauginos producing 
the trilepton signal have masses of m-± ~ m^o = 184 GeV, and the heavier gauginos, which 
decay mostly into the lighter gauginos, W and Z bosons as well as the lightest Higgs boson, 
have masses of m^o = 400 GeV and m-± ~ m^o = 415 GeV. The average squark and gluino 
masses are rriq ~ 550 GeV and m~ g = 604 GeV. Note that this benchmark point is also 
relatively close to the region excluded by the Tevatron collaborations CDF and DO, which 
assume, however, a lower value of tan/3 = 3 and A = {2]. 



A. Transverse momentum spectra 

In Fig. [2J we present transverse momentum spectra for light charged and neutral gaugino 
pairs with masses of 113 and 184 GeV at the Tevatron (top) and early LHC (bottom), where 
the center-of-mass energies are y/S = 1.96 and 7 TeV, respectively. We show predictions 
at fixed order 0{oc s ) (dashed) as well as with transverse-momentum (dotted) and joint 
resummation (full curves). While the fixed-order predictions diverge at small px due to an 
uncancelled soft singularity from real gluon emission, the resummed predictions exhibit a 
finite, physical behavior with a pronounced maximum in the region of pt = 5 to 10 GeV. In 
the region of intermediate pt of 20 to 60 GeV, the resummed predictions are considerably 



11 



p p ^ % % at the Tevatron (1.96 TeV) 

i j 




20 40 60 80 100 

P T [ GeV ] 



P P -» XX at the LHC (7 TeV) 




20 40 60 80 100 

P T [ GeV ] 

FIG. 2: Transverse momentum distributions of light gaugino pairs at the Tevatron (top) and LHC 
(bottom) with = 1.96 and 7 TeV center-of-mass energy, respectively, in fixed order (dashed) 
as well as with transverse-momentum (dotted) and joint (full curves) resummation. 
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larger than those at fixed order. The calculations using transverse-momentum and joint 
resummation are in good agreement, but the theoretical uncertainty of the latter, estimated 
by varying the renormalization and factorization scales by a factor of two about the average 
mass of the two gauginos, is considerably smaller, since threshold logarithms are resummed 
simultaneously. The cross sections for chargino pairs exceed those for the trilepton channel 
at the Tevatron and for the X2X1 channel (but not the X2X1 channel, which is not shown) 
at the early LHC with the cross section for pair production of the second-lightest neutralino 
being more than an order of magnitude smaller at the early LHC. 

At the LHC design luminosity of yS = 14 TeV, the neutralino pair production cross 
section is larger by more than a factor of three (see the top part of Fig. [3]). Otherwise, the 
behavior of the various predictions is very similar to the one described above. One notices, 
however, that the resummed predictions exceed those at fixed order only at larger values of 
Pt > 40 GeV. Furthermore, the trilepton cross section for positive charginos is larger than 
the one for negative charginos by almost a factor of two, since in contrast to the Tevatron 
the LHC is a pp collider. With a center-of-mass energy of VS = 14 TeV, it may also become 
possible to observe the associated production of heavier neutralinos xl an d charginos xf 
with masses of about 415 GeV. The corresponding transverse momentum spectra are shown 
in the lower part of Fig. [31 The absolute values of the cross sections are reduced by about 
a factor of 50 for both positive and negative charginos, but the shape of the distributions is 
very similar. Joint resummation leads again to the smallest scale uncertainties. 

B. Invariant mass spectra 

While the joint resummation formalism is designed to match more closely the one for 
transverse-momentum resummation, it also allows to simultaneously resum threshold log- 
arithms and obtain precise invariant mass spectra. These are therefore presented in this 
section for various gaugino pairs and colliders and compared to those obtained with pure 
threshold resummation with the expectation that the agreement will be slightly worse than 
the one for transverse momentum distributions. 

In Fig. H] we show invariant mass spectra for the trilepton channel at the Tevatron obtained 
at LO (short-dashed) and with NLO SUSY-QCD corrections (dashed) as well as with pure 
threshold (dotted) and joint (full curve) resummation. In this figure, the NLO and threshold 
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P P -> XX at the LHC (14 TeV) 

i j 




20 40 60 80 100 

P T [ GeV ] 

P P -+% 4 jt 2 at the LHC (14 TeV) 

Joint resummation 



P T resummation 




20 40 60 80 100 

P T [ GeV ] 

FIG. 3: Transverse momentum distributions of light (top) and heavy (bottom) gaugino pairs at the 
LHC with y/S = 14 TeV center-of-mass energy in fixed order (dashed) as well as with transverse- 
momentum (dotted) and joint (full curves) resummation. 
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M [ GeV ] 

FIG. 4: Invariant mass distributions of light gaugino pairs at the Tevatron with y/S = 1.96 TeV 
center-of-mass energy in fixed order (short-dashed, dashed) as well as with threshold (dotted) and 
joint (full curve) resummation. 

resummed predictions are both considerably larger than the one obtained at LO. In the linear 
representation of da/dM emphasizing the low- invariant mass region shown here, they can in 
fact not be distinguished, as threshold effects only start to dominate as the invariant mass 
squared M 2 approaches the total available center-of-mass energy -^/s. This is also the reason 
why the jointly resummed prediction differs and in fact exceeds slightly the pure threshold 
resummed prediction at small M, as large logarithms at small have been simultaneously 
resummed. This leads also to an additional reduction of the scale uncertainty, represented 
again as a shaded band and obtained by varying the renormalization and factorization scales 
by a factor of two about the invariant mass M. 

The various features described above are even more prominent at the LHC with its larger 
design center-of-mass energy of v^S* = 14 TeV, despite the fact that also the masses of the 
light gauginos are slightly larger in the SPSla' scenario than at the LMO benchmark point. 
The scale uncertainties in the upper part of Fig. [5] are considerably smaller at LO and even 
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FIG. 5: Invariant mass distributions of light (top) and heavy (bottom) gaugino pairs at the LHC 
with y/S = 14 TeV center-of-mass energy in fixed order (short-dashed, dashed) as well as with 
threshold (dotted) and joint (full) resummation. 
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more at NLO and NLL+NLO, as the variation from M/2 to 2M is less important compared 
to the large yS. This leads to threshold and joint resummation predictions that no longer 
overlap within scale uncertainties in the region of maximal cross section, indicating that 
the theoretical error is slightly underestimated in this case. As it was already observed in 
the previous section, the cross section for the trilepton channel with positive charge exceeds 
the negative one by almost a factor of two due to the positively charged initial state at the 
LHC. Heavier gauginos are produced closer to threshold, so that the distributions in the 
lower part of Fig. [5] are found at larger values of M. The scale variation and the shaded 
bands representing it become more important again with respect to \/S as do the threshold 
logarithms, so that the threshold and jointly resummed predictions overlap again within the 
theoretical uncertainties. The absolute size of the cross section is again almost two times 
larger for X4X2 pairs than for X4X2 pairs. 

C. Total cross sections 

Total cross sections can be obtained from the distributions shown in the previous sections 
by integrating over either p? or M. This is numerically not always trivial, but should in 
principle lead to similar predictions, building confidence in the theoretical calculations. In 
particular, since the perturbative, resummed and expanded contributions are obtained in 
(Pt, M) and (b, N) space, respectively, we first compute the differential cross section da/dpr 
and then integrate over p? with the trapezoidal rule only after the different contributions 
to the total cross section have been matched. 

The results are displayed in Tab. [H They show the expected important increase in 
absolute size and reduction in scale uncertainty from LO to NLO and then, to a lesser 
extent, at NLL+NLO. It is interesting to note that the pr-resummed predictions are indeed 
larger than those obtained at NLO, but the threshold resummed predictions are slightly 
smaller. The jointly resummed predictions are very similar to those obtained at NLO at the 
two lower center-of-mass energies and larger at the LHC with \/S = 14 TeV. As expected, the 
lowest scale uncertainties are found with threshold resummation when the particle masses 
are large compared to the available center-of-mass energy (light gauginos at the Tevatron, 
heavy gauginos at the LHC with v^S* = 14 TeV), while joint resummation gives results that 
are similar to pt resummation, but more precise, in the two other cases. 
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TABLE I: Total cross sections (in fb) with scale uncertainties for the trilepton channel at the 
Tevatron in the LMO scenario and for light and heavy gaugino pairs at the LHC with 7 and 14 
TeV center-of-mass energy in the SPS la' scenario at various levels of accuracy. 



Collider 



Gauginos 



LO 



NLO 



PT 



Threshold 



Joint 



Tevatron 
LHC-7 
LHC-14 
LHC-14 



X2xf 
x\X\ 
X2X1 
X4X2 



300.0441^-555 3 66 .794+io.772 3 7 7. 118 +5.89i g 363.922+^02 365.974+3.542 
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121.216; 



-1.843 
-1.557 



122.81711;^ 
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346.5381^1° 419.930™ 429.678l°;j*§ 416.327l°|^ 428.2021°;** 
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7.844l° : ^ 



8.261±g:i™ 



7 7«Q+0.017 
' • ' uo -0.008 



o 1 1 7+0.033 
°- LL ' -0.047 



IV. CONCLUSION 



In this paper, we have completed our investigation of gaugino production at hadron 
colliders with different resummation methods by presenting a NLL+NLO calculation that 
jointly resums large logarithms in the small-pr and threshold regions. After a detailed 
outline of the organization of the analytical calculation and the numerical implementation, 
in particular of the convolutions with the parton densities and the inverse Fourier and Mellin 
transforms, we have compared the new jointly resummed predictions to those obtained 
previously at fixed order as well as with pure px and threshold resummation. We found 
in general good agreement in the transverse momentum and invariant mass distributions, 
confirming the reliability of the joint resummation method. In most cases the new results 
were also even more precise, so that they can be considered the most reliable predictions 
for direct gaugino production available for the Tevatron and the LHC. They are therefore of 
great importance for supersymmetry searches and parameter determinations at these hadron 
colliders and should be taken as the basis for future experimental analyses. 
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